Radiation Damping in FRW Space-times with Different Topologies 
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We study the role played by the compactness and the degree of connectedness in the time evolution 
of the energy of a radiating system in the Friedmann-Robertson- Walker (FRW) space-times whose 
t — const spacelike sections are the Euclidean 3-manifold 1Z 3 and six topologically non-equivalent 
flat orientable compact multiply connected Riemannian 3-manifolds. An exponential damping of the 
energy E(t) is present in the 1Z 3 case, whereas for the six compact flat 3-spaces it is found basically 
the same pattern for the evolution of the energy, namely relative minima and maxima occurring at 
different times (depending on the degree of connectedness) followed by a growth of E(t). Likely 
reasons for this divergent behavior of E(t) in these compact flat 3-manifolds are discussed and 
further developments are indicated. A misinterpretation of Wolf's results regarding one of the six 
orientable compact flat 3-manifolds is also indicated and rectified. 



PACS numbers: 98.80.-k, 98.80.Hw, 04.20. Gz, 04.20.Jb, 04.20.Cv 
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As general relativity is a purely metrical (local) theory it clearly leaves unsettled the global structure (topology) 
of space-time. However, in cosmology perhaps the most important problems are related to the global structure of 
space-time, where the topological degrees of freedom ought to play an essential role. 

Geometry constrains, but does not determine the topology of a space-time. Consider, for example, the Friedmann- 
Robertson- Walker (FRW) space-times, whose line element can be given by 



ds 2 = dt 2 - A 2 {t) 



dr 2 



+ r 2 (d9 2 + sin 2 9 dip 2 ) 



(1.1) 



where A(t) is the scale factor, t is the cosmic time, and the constant spatial curvature k = 0, ±1 specifies the type of 
geometry (flat, elliptic or hyperbolic) of the t — const spacelike section M3. Clearly FRW space-time manifolds M4 
can be splitted into KxM 3 . The number of three-dimensional spacelike manifolds M3 which can be endowed with 
the three possible geometries of the t = const 3-spaces of FRW space-time is quite large jl],@]: for k = there are 18 
topologically distinct 3-spaces, while for both k = ±1 an infinite number of 3-spaces exist ^,0]. Even if we restrict 
ourselves to orientable and compact manifolds we still have an infinite number of 3-spaces At 3 for the elliptic and 
hyperbolic cases, and six families of topologically different spacelike manifolds for the flat 3-spaces - . 

Since physical laws are usually expressed in terms of local differential equations, in order to be confident about 
the physical results one derives it is often necessary to have some degree of control over the topological structure of 
the space-time manifold so as to include constraints imposed by the topology @]. One is then confronted with the 
question of what topologies are physically acceptable for a given space-time geometry. An approach to this problem is 
to study the possible observational (or physical) consequences of adopting particular topologies for the space-time 0] 
- Jr| (see also || and references therein). 

In this work we study the role played by the compactness and the degree of connectedness in the time evolution of 
the energy of a radiating system in the flat FRW space-times whose spacelike t = const sections are endowed with 
seven different topologies, namely the simply connected 3-space 1Z 3 , and six multiply connected orientable compact 
3-manifolds obtained by suitable identifications of opposite faces of cubes (four) and of hexagonal prisms (two) after 
suitable turns ||,|,|lj- 

The radiating system we shall be concerned with is represented by a point-like harmonic oscillator (energy source) 
coupled with a relativistic massless scalar field |2C]] - (^9|. Similar radiating systems have been used to study a 
wide class of radiation phenomena as, for example, gravitational waves [ p3| or the radiated energy of oscillating 
electromagnetic dipoles 30 31j] . Our radiating system was also used in the study of dynamic and geometric constraints 
on the radiation in elliptic FRW expanding universes [E7J - M . 
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In the next section the radiating system is described: the action integral is presented, its variation is performed and 
the cor resp onding evolution equations are obtained. We also outline there a proposal for solving these equations. In 
section [ill] the Green functions of the wave operator are derived for each specific space-time manifold with arbitrary 
scale factor A(t), and combined with the evolution equations to give the corresponding radiation reaction equations. 
A misinterpretation M - 0] of Wolf's Q results (Theorem 3.5.5) regarding one of the six orientable compact flat 



3-manifolds is also indicated and rectified therein. In section IV, the radiation reaction equations are numerically 
integrated. Graphs that show how the system energy varies with the time for each topologically different flat FRW 
manifold are presented for static and non-static cases. Our main conclusions are discussed in section |v|. We show 
that when the t = const sections are the simply connected 3-space 1Z 3 the radiation damping phenomenon is present, 
whereas for all compact flat FRW space-time manifolds we have investigated, the energy E(t) exhibits a few relative 
minima and maxima followed by a growth of the energy with the time. Possible reasons for this divergent behavior of 
E{t) in these compact flat 3-manifolds are examined. We also discuss the role played by the compactness and by the 
degree of connectedness in the energy patterns of our system in these space-time manifolds. Further developments 
are indicated therein. 



II. PHYSICAL SYSTEM AND EVOLUTION EQUATIONS 



Radiation waves produced by an oscillating energy source as, e.g., an oscillating electromagnetic dipole, have been 
studied by using a theoretical model represented by a classical harmonic oscillator (e nerg y source) coupled with a 
relativistic massless scalar field (scalar radiation waves propagating at speed of light) [[20| - [ p8| . In our model the 
gravitational field is treated as external, but a suitable conformal coupling with the scalar field is considered ]32] ]. 

The dynamics of our system can be described by an action integral, which contains a term for the scalar field 
<p(t, x), another associated to the oscillation amplitude Q(t) of the point-like harmonic oscillator, and a coupling term 
between the scalar field and the harmonic oscillator according to 

S=\J d i x^[g^d^d u( j>-^R ( t> 2 ] + \ J dt[Q 2 -cj 2 Q 2 } 
+ A / d 4 x p(t, x) Q(t) cf){t, x) , 



(2.1) 

where t € [to, oo), x £ M3, g^ v is the metric tensor on VW4, g = det [g^v) , R is the scalar curvature of VW4, overdot 
means derivative with respect to t, and A is a coupling constant. The function p is the normalized density function, 
which accounts for the coupling between the harmonic oscillator and the scalar field. Similarly to the coupling between 
charges and electromagnetic fields in classical electrodynamics, we shall consider in this work a point-like coupling 
between the harmonic oscillator and the scalar field, namely the one in which p(t,x) = 5^ (x) / \/—g(t, x), where 
is the three-dimensional Dirac delta function. This type of coupling requires a renormalization of the frequencies, and 
to this end we need an e-family of uncoupled frequencies lo £ (see section III). Here and in what follows units in which 
c = 1 are used. 

Varying the action (2.1) with respect to 4> and Q one obtains the coupled evolution equations of the system, namely 







<p(t,x) = Xp(t,x)Q(t), 



Q(t)+uj 2 (t)Q(t) = A / d?x«J=g~p{t,g)<t>(t,x), 



(2.2) 
(2.3) 

where □</> = {*/—■ ff) 1 d^^J—g g^ v d v cf)) is the d'Alembertian operator, and □ + ^R is the wave operator defined on 
Ma, hereafter simply called wave operator |P2| , p3[ . 

Before proceeding to the discussion of the Green functions for the wave operator in flat FRW space-time mani- 
folds (for any smoot h A (i)) we shall consider how o ne can obtain the radiation reaction equation from the above 
equations (2.2) and ( [2.3] ). One first solves equation ( |2.2| ) as an initial value problem, by writing the solution in the 
form H|] <fr(t,x) — (j>i(t,x) + <pH{t,x), where (j)H(t,x) satisfies the corresponding homogeneous equation, and where 
the solution of the inhomogeneous equation is given by 



fa(t,x)= / dt' d 3 x' y/-g(t', x') G(t, x; t' , x) A p(t', x) Q(t') , 



with t' G [to,*], x' € Mz- In 
solution of the wave operator |33[ 



(2.4) 



.41) Q(t, x;t' , x 1 ) is the retarded Green function, often referred to as fundamental 
Since the scalar field of our system may be thought of as the propagation medium 
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for the radiating energy, we assume the initial condition (<f)(to,x),dt<f>(to,x)) = (0,0), which means that the scalar 
field carries no energy at t = to. This condition implies that 4>ij(t, x) = 0, and therefore (j>(t, x) = 4>i(t, x). Finall y, th e 
radiation reaction equation of the energy source Q{t) can be found by using the relation <fi = <j> [Q] in equation fl2.5p. 
We shall return to this point later in the next section. 

In this work, the Green function plays an essential role in that it incorporates both the geometrical and topological 
features of the FRW t = const 3-spaces. They will be obtained, in the next section, through the study of null geodesies 
of the space-times for each distinct 3-space (see table Q). 

III. GREEN FUNCTIONS FOR THE WAVE OPERATOR 

In this section we shall derive the Green functions of the wave operator for the flat FRW space-times whose spacelike 
t = const sections are the multiply connected compact orientable 3-spaces, obtained by suitable identification of faces 
of a basic cell (see Refs. [HJJ]) according to the table |, and the simply connected space 1Z 3 . Clearly the metric tensor 
we shall be concerned is given in cartesian coordinates by 

9tlv = diag ( 1, -A 2 (t), -A 2 (t), ~A 2 (t) ) . (3.1) 

Regarding the manifold one often encounters, in the literature on compact orientable flat 3-spaces, reference to 
a cube in which each pair of opposite faces is identified after half a turn - Hj. However, if such a cube is endowed 
with the Euclidean geometry then the resulting three-space is not a manifold, but an orbifold p^ , ^6| , ^i| . This 3-space 
fails to be a manifold along the edges of the cube. Only when the cube is endowed with the elliptic geometry the 
resulting space is indeed a manifold, namely the well known real projective space V 3 . 

Recently, using theorem 3.5.5 by Wolf JO], it has been shown that a basic cell for the Euclidean manifold T4 is the 
cube shown in figure 1 (4^] . It should be stressed that the basic cell for T4 shown in figure 1 is a standard basic cell; 
in general, however, the basic cell for T4 needs not to be a cube. 



Topology Type 


A Basic Cell 


Identifications of Faces 


Ti 


cube 


3 pairs non rotated 


T 2 


cube 


2 pairs non rotated, 1 pair rotated 90° 


% 


cube 


2 pairs non rotated, 1 pair rotated 180° 


T 4 


cube 


1 pair rotated 180°, 2 pairs according to figure 1 


Hi 


hexagonal prism 


top and bottom rotated 60° 


H2 


hexagonal prism 


top and bottom rotated 120° 



TABLE I. The six compact orientable topologies for the flat 3-manifolds can be obtained by identifying opposite faces of a 
basic cell as shown in this table. 
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The cubes and TZ 3 cases 

We shall discuss in this part the Green functions for the wave operator on 3-spaces % (i = 1, • • • ,4) shown in 
table |. The manifold 71, referred to in the literature as the three-torus T 3 p7§ , is a compact multiply connected 
Riemannian manifold, which can be obtained by identifying the opposite faces of a cube of side a. 

To build the Green function in TZ x 7i we first use the conformal properties of FRW space-times. Defining the 
conformal time r = f(t) = Jdt'/A(t') for t 1 6 [to,t], we transform the problem of finding the Green function 
G(t, x;t' , x 1 ) of the operator (□ + \K) m to the problem of finding the Green function G(t,x;t' = 0,2?') of the 
operator □ |^2|,^3|. Moreover, a useful way of thinking about 71 in terms of the simply connected 3-manifold 1Z 3 is 
to imagine the cube repeated endlessly in a three-dimensional grid (basic cell and its images), where each repetition 
consists of the same physical region of space. The space TZ 3 containing an infinite grid of copies of the basic cell is 
called the universal covering space of l~i . In this way we now transform the problem of calculating the Green function 
of □ at (r, x) € TZ X T\ due to a point-like source located at x! = (0,0,0) and emitting at r' = 0, to the equivalent 
problem of finding the field at (r, x) € TZ x TZ 3 due to an infinite set of point-like sources, each one located at the 
center of a cube of the grid, and all of them irradiating simultaneously at r' = 0. 

In this infinite grid picture the Green function of □ is calculated summing up the contributions of each one of the 
point-like sources. Since a scalar ray is a null geodesic, the distance travelled by each scalar wave front is a measure 
of the corresponding time travel, so the conformal time r of arrival at x will depend on the position at which each 
source is located. Null geodesies corresponding to retarded waves connecting (r' = 0, a(n)) to the observation point 
(t, x) are then such that r — \x — a(n)\ = 0, so the Green function for r > is given by 



4ttR(t)R(t' = 0) \x-a(n)\ 



where x, a(n) S TZ 3 , a(n) = a(n\ei + n2&2 + n^e^) . Here {e{\ is the usual orthonormal basis vectors in TZ 3 and 
(711,712,713) are integers, and R(t) = A(/ _1 (r)). 
We s hall consider now the radiation reaction eq uati on of our energy source in T\, Using the Green function 



G and fl2.4| ) one finds ef>[Q], which can be used in (2__3) to furnish the radiation reaction equation of our system. 
Note that in these calculations the first term (77,1 = 712 = 13 = 0) of the Green function gives rise to the term 
(A 2 /4 7r) J dxS(x)/(A(t) x), that formally diverges. A suitable renormalization procedure is therefore necessary. This 
resembles the need for renormalization one finds when dealing with accelerated point charges in classical electrody- 
namics [ pCfl . We learn from p8| that a renormalization procedure can always be made in this case. This amounts to 
saying that for any < Q 2 < 00, and any e > 0, one can always define an e-parameter family of frequencies ui 2 by 

^Etftl/^, (3.3) 

such that p e (x) is a ^-sequence that in the limit e — > converges to the Dirac delta 8(x). 

Using now (|3.3| ), in which we have made 2T = A 2 /(Air) , one finds that, for arbitrary initial data (Q(to), Q(to)), the 
renormalized radiation reaction equation of the source can be written in the form 

Q(t) + 2TQ(t) + n 2 Q(t) = ^- (3-4) 

^ ' Til, "2, "3 

for any t € [to, 00), and where a(n) = = a [n\ 2 + 712 2 + 713 2 ] 1 / 2 for all (711,712,713) integers not all zero, 

t n = f^ 1 (a(n)) and the step distribution O(s) = for any s < 0, O(s) = 1 for all s > 0. Note that the radiation 



reaction equation (3.4) contains an infinite, but countable, number of retarded terms. Note also that, for a given pair 



of initial data (Q(to), Q(to) ), the continuity of Q and Q for all t is sufficient to ensure that (3^4) can be integrated to 
give a unique solution. 

The next 3- manifolds we shall consider are T2, T% and T± (see table |). The 3-space T2 is a compact multiply 
connected flat three-dimensional Riemannian manifold, obtained by identifying opposites faces of a cube, with a pair 
being identified after a rotation of 90° of one face relative to the opposite face. The 3-manifold T% is obtained by 
identification of opposites faces of a cubic cell, but now a pair of faces are identified after a turn of 180° (see table |). 
Finally, the 3-manifold is obtained by pairwise identifying the faces of the cube according to figure 1. For these 
3-spaces the Green functions for the wave operator can again be built by using the infinite grid of cubes picture. It 
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should be noticed, however, that the grids of cubes corresponding to T%, T 2 , T 3 and T4 are different []. Nevertheless, 
as the center of a basic cell in each case represents the same infinite set of points (images of the point-like source at 
(0, 0, 0) ) of the covering space 1Z 3 , the Green functions for the wave operator in these 3-spaces are equal to the one 



we have calculated for T\, i.e., are giv en b y (3.2). Thus, for arbitrary initial data (Q(to), Q{to)) the radiation reaction 
equation in T2, T3 and T4 reduces to (3.4) obtained for 7±. 

A word of clarification is in order here: the above Green functions for the wave operator for the 3-spaces % 
(i = 1, • • • , 4) coincide only because in each case the point-like source is located at the center of the basic cell x' = 0. 
This is so because although all these 3-spaces are locally homogeneous, the 3-spaces T2, T3 and T4 are not globally 
homogeneous, thus the corresponding Green functions of □ for x 1 7^ are more involved and do not depend simply 
on the relative position between the source and the observer. 

Regarding the simply connected flat 3-dimensional Riemannian manifold TZ 3 (t = const section of the Minkowski 
space-time), the Green function for the wave operator □ is equal to the first term (n\ = 112 — H3 = 0) of G appearing 



in ( J3 - 2| ) - Thus the radiation reaction equation in this manifold reduces to (3.4) with right hand side equal to zero, for 
all t G [t , 00). 

The hexagonal prism cases 

We shall discuss now the Green functions for the wave operator on the flat hexagonal prism manifolds Ti.\ and TI2 
described in table Q. 

The 3-space Hi is a compact multiply connected flat three-dimensional Riemannian manifold, which can be obtained 
by identifying the top and bottom faces of a regular hexagonal prism after a rotation of 60°, while the lateral faces 
are pairwise identified in the usual manner. The 3-manifold H2 can be similarly defined, but now the identification of 
the top and bottom faces (regular hexagons) is made after a turn of 120°. In both cases we denote by a the shortest 
distance between two opposite sides of the regular hexagon and by h the height of the hexagonal prism cell. 

For a point-like source located at the center of the basic cell, the Green functions for the wave operator in the 
3- manifolds Hi and H2 are obviously the same. Moreover they turn out to be similar to the Green function (ft.2f). 
But now instead of a(n) one has b(n) = a n\^-e\ + a {\nx + U2)&2 + h n^e^ , where (m , n.3) are integers, to locate 
the centers of the hexagonal prism cells in which TZ 3 has been tessellated. Note that the number of images at a given 
distance from the point-like source clearly depends upon the ratio h/a. The radiation reaction equation in the present 



cases is similar to (3.4) but with the obvious change a(n) — > b(n) = \b(n)\ = [a 2 (ni 2 + n\7i2 + ri2 2 ) + h 2 ^ 2 } 1 / 2 and 
again the integers (ni, 712,713) are not all zero. 

It should be emphasized that as the 3-manifolds Tlx and H2 are locally but not globally homogeneous, the Green 



function for the wave operator in these manifolds is again given by (3^2), with b(n) instead of a(n), only because the 
point-like source is at the center of the basic cell x 1 = 0. 

To close this section we emphasize that the Green functions, obtained for flat FRW space-time manifolds M4 = 
TZ x contain the topological constraints of .M3, information that makes possible to find out the exact radiation 
reaction equation in each case. 

IV. NUMERICAL ANALYSIS 

In this section we shall discuss the time behavior of the energy of the harmonic oscillator E{t) = \ [Q 2 (t) + ft 2 Q 2 (t)] , 
where the function Q(t) is the solution of the radiation reaction equation corresponding to each flat manifold we have 
discussed in the previous section. 

Without loss of generality, in the integration of the radiation reaction equations and in the plotting of the energy 
function we have taken specific values for the constants. We have also assumed the continuity of Q and Q, and 
chosen suitable values for the initial data Q{to) and Q(to). For a neat comparison between the simply and multiply 
connected cases we have chosen r = 1, ft 2 = 30, the length a = 1, and for the heights of the hexagonal prism and 
three-torus we have taken h — 0.4 and 1. As a matter of fact, our three-torus T\ w as obtained from a parallelepiped 



with edges a, a, h, so a(n) — [a 2 (ni 2 + n,2 2 ) + /i 2 ^ 2 ] 1 / 2 was used in equation (3.4). Further, we have also chosen as 
initial data t = 0, (Q(0),Q(0)) = (^2/0,0), which means that the initial energy of the source has been normalized, 
i.e., 75(0) = 1. 



1 This is so because each 3-manifold is obtained by forming the quotient Mg — 1Z 3 /Fi (i = 1, ■ • ■ , 4), where for each case T; 
a different discrete group of isometries of the covering space 1Z 3 acting properly discontinuously, without fixed points |H|. 
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Figures 2, 3 and 4 correspond to the static case A(t) = const = Ao, which we have normalized to 1 (Aq = 1), while 
dynamic situations (A(t) ^ 0) are considered in figure 5. 

Taking into account the above choices of values and using the computer algebra systems Mathematica |38| and 
Maple |39| ] the numerical integrations of the radiation reaction equations as well as the corresponding graphs for the 
energy function were obtained (see figures 2 to 5). 

Figure 2 shows the behavior of the energy with the time for the flat FRW space-times, in which the t — const 
sections are any of the orientable compact multiply connected 3-manifolds 71, 7~2 and T3 and for different ratios 
h/a. Note that h is the height of the parallelepiped of square basis of side a. The curves exhibit basically the 
same pattern, namely relative minima and maxima followed by a predominant growth of the energy with the time 
(E(t) — > 00 when t — ► 00). These extrema are related to the contribution of the discrete retarded terms (the right hand 
side of the radiation reaction equation) demanded by the compactness and the corresponding connectedness of these 
manifolds. The fact that the relative extrema occur at different times and are of different amplitudes (intensities) for 
distinct tessellations (different ratios h/a) of the covering manifold, basically reveals the differences in their degree of 
connectedness (returning rays take different times to return to the origin). This growth of the energy is discussed in 
the next section. 

Regarding the behavior of the energy function for the cases in which the t = const sections are Hi and H2, shown 
in figure 3, we again note that due to the compactness and connectedness we have relative minima and maxima with 
distinct intensities, which take place at different instants for distinct ratios h/a. The curves again display an eventual 
growth of the energy with the time for these manifolds. 

Figure 4 compares the variation of the energy E with the time for the cases where the t = const section M 3 are 
7Z 3 , T[ (with a = h = 1) and Hi (also with a = h = 1). This figure shows for the Minkowski space-time, as expected, 
an exponential decay of the energy with the time, whereas for the manifolds T\ and Hi it shows basically the same 
pattern, i.e. relative minima and maxima occurring at different times, depending on the degree of connectedness, 
followed by a growth of the energy with the time. 

Although the net role played by the degree of connectedness as well as compactness can be singled out in the static 
cases A(t) = const = A$, for the sake of completeness we have examined three instances where dynamic expansion 
takes place. Figure 5 corresponds to the plot of the energy function for the Ti manifold (with a = h = 1) in the 
dynamic expanding cases: (i) linear expansion A(t) = t + 0.7, (ii) square root expansion A(t) — \ft 1 (iii) inflationary 
expansion A(t) = e*/^, that is, an expansion with future event horizon. 

Although for the case (iii) one clearly has radiation damping, we have not been able to find out so far a closed 
formal proof of the asymptotical behavior for the other two cases. We emphasize, nevertheless, that the net role played 
by the connectedness and compactness can be better singled out in the static cases, where the dynamical degrees of 
freedom are frozen. The study of that role for the static cases is in fact the major aim of the present work. 



V. CONCLUSIONS AND FINAL REMARKS 



In this work we have studied the role played by the topological compactness and connectedness in the time evolution 
of the energy of an harmonic oscillator in flat FRW space-time manifolds, whose t = const sections are (i) the orientable 
simply connected non-compact 3-space TZ 3 and (ii) six possible flat orientable multiply connected compact 3-manifolds 
given in table |. 

For the 1Z 3 case we found that the energy function E(t) exhibits an exponential decay with the time — the radiation 
damping (E(t) — ► when t — ► 00) takes place, as one could have expected in agreement with |2(J - p6[ . 

For the manifolds 71, T2, T3 and T4 as well as for the manifolds Hi and H2 the behavior of the energy with the 
time exhibits basically the same pattern: relative minima and maxima occur at different times for distinct ratios h/a 
(distinct tessellations) depending on the degree of connectedness of each 3-manifold, and are followed by a growth of 

This asymptotical divergent behavior of E(t) for these compact manifolds contrasts with the radiation damping of 
the energy we have found for 1Z 3 . There is a quite simple heuristic argument which suports our numerical results, 
though. If one ad hoc assumes an exponential asymptotical behavior for Q, i.e. Q(t) = 7 exp ((3t) with (3 and 7 real 
constants, then for the static cases (A(t) = 1), and for each of the above compact manifolds, in the limit t — > 00 
equation ( |Q| ) reduces to 

/? 2 + 2r/? + r> 2 ~ Cm 

= }^ — exp(-/3a m ) , (5.1) 

z 1 a m 

m—l 

where c m is the number of images of the point-like source at a distance a m in the infinite grid picture. To attain our 
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goal, w e wi ll show that (5.1) has only one real solution for f3, which is positive. Indeed, let f((3) be the right hand 



side of (5T), which is a positive monotone decreasing function of /3, and such that 



lim /(/?) = oo and lim /(/?) = . (5.2) 

Thus f((3) lies entirely in the first quadrant of the plane a nd c rosses it from the top-left to the bottom-right. Now, 
since T > then for a given pair (r, fl) the left hand side of (|5.l|) is a parabola curved upwards with vertex at f3 = —V. 
Therefore, it always intersects the curve for f([3 ) in just one point, which is in the first quadrant. In other words, 



there is only one real (3 solution to equation (5.1), which is positive. 

The unexpected (unphysical ?) growth of the energy with the time for the above compact flat 3-manifolds cases 
illustrates that non-trivial topologies can induce rather important dynamic changes in the behavior of a physical 
system. This type of sensitivity has been refered to as topological fragility and can occur without violation of any 
local physical law |lq |. A rigorous non-numerical analysis of the reasons for this surprising divergent behavior of E(t) 
when compact flat FRW space-times are considered has been carried out, and we hope to publish our results shortly 
elsewhere. We anticipate, however, that the causes for such a behavior lie in the compactness of the manifold in at 
least one direction, on the one hand, and in the type of coupling between Q and <f>, on the other hand. 

A possible physical measure of the degree of connectedness in these compact 3-manifolds can be made through the 
study of the number of emitted rays that return to the origin within a given lapse of time. According to this concept 
of degree of connectedness one learns from figures 2, 3 and 4 that the greater is the degree of connectedness the earlier 
is the occurrence of the first relative minimum in the energy function. Incidentally, note that the extension of this 
concept of degree of connectedness to non-compact 3-manifolds implies that TZ 3 has a null degree of connectedness. 
This, of course, is indicated in figure 4, which shows a net exponential decay of the energy with the time; no relative 
minima and maxima come about, which means that no ray returns to the origin. 

A simple inspection of the graphs for E(t) clearly shows that the derivative E — Q(0 + il 2 Q ) of the energy 



function is discontinuous at a few points. Indeed, for the static case, for example, using (3.4) one obtains 

E{t) = 2YQ{t) V — T Q(t - a(n)) Q(t - a{H) ) - Q(t) 
^-^ a(n) 

n\,n-2,n 3 x ' 



(5.3) 



From this equation one sees that the discontinuities occur at t = a(n) , th at is, they come about each time a new 
term Q(t — a(n)) / a(n) is taken into account in the right hand side of (|3.4[) . A question which naturally arises here 
is whether the inverse problem, i.e. that of determining the basic cell (topology) corresponding to the spacelike 
t = const sections from the graphs of the energy E(t), can be solved. Regarding this problem it is clear that one can 
find the distances of the point-like source to its images by using where the discontinuities of E(t) take place, and the 
number of images at a given distance through the magnitude of the corresponding discontinuities. So, one can probe 
the topology of the 3-spaces at least in a few cases. It is not yet clear whether an algorithm for solving the inverse 
problem for the most general (Euclidean) setting can be found, though. 

As far as we are aware || this is the first work in which a physical consequence of adopting the flat hexagonal 
prisms Tii and TL2 has been studied. 

It is of worth to emphasize that when the expansion of the universe is considered the degree of connectedness is less 
than the ones for the static cases. For the manifold T\ this can be seen by comparing figure 2 (static case A(t) = 0) 
and figure 5 (monotone expansions A(t) > 0). 

Before closing this article we should like to stress that our study does not cover all possible spatially compact 
oricntable flat FRW manifolds. Thus, for example, we have not considered that for the 3-manifolds T\ and T% the 
basis of the basic cell need not to be a square, it can be a parallelogram. The restriction we have made, however, does 
not seem to be decisive for the patterns of the behavior of the energy with the time we have found. 

To conclude we remark that the study of radiation damping in elliptic (k = 1) FRW manifolds in which the t = const 
sections are endowed with different (orientable compact) topologies is being carried out. 
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FIG. 1. A basic cell for the manifold Ta is a cube whose faces are pairwise identified as indicated in this figure. 

FIG. 2. Behavior of the energy of the harmonic oscillator for T = 1 and Q 2 — 30 in static flat FRW space-times with T\ 
space slices. There are a few relative maxima followed by a growth of E(i). Two different ratios h/a are considered (h is the 
height and a is the side of the square basis of the basic cell). 
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FIG. 3. Behavior of the energy E(t) of the harmonic oscillator for T — 1 and Q 2 = 30 in static flat FRW space-times with 
3-space Hi. There are a few relative maxima followed by a growth of E(t). It shows the energy vs. time curves for distinct 
ratios h/a (h is the height of the hexagonal prism and a is shortest distance between two opposite sides of the regular hexagon). 

FIG. 4. The time evolution of the energy of the harmonic oscillator for flat, static FRW space-times with different topologies 
for the spacelike t = const sections: TZ 3 , T\ and H\ (both with a — h = 1). Different degree of connectedness implies different 
patterns of E(t). Here again F = 1 and Q 2 — 30. 

FIG. 5. The time evolution of the energy of the harmonic oscillator in FRW expanding space-times with 3-space 71 (for 
a = h = 1). Three types of dynamic expansion are shown: linear [A(t) = 1], square root [A(t) = t 1 ^ 2 ] and inflationary 
[A(t) = ae at , a' 1 = 2 1/2 ]. Here again T = 1 and ft 2 = 30. 
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